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(54) Fat free NMR TOF angiography 

(57) A method for performing motion -compensated 
spectral-spatial selective magnetic resonance imaging. 
The method includes applying a radio frequency (RF) 
spectral-spatial excitation pulse sequence to a region of 
a body. The excitation sequence includes a binomial 
pulse train including at least two sub-pulses. The bino- 
mial pulse train also provides spectrally selective exci- 
tation of a selected resonance frequency at the prede- 
termined field strength. The method also includes ap- 
plying to the same body region an oscillating slice se- 
lection gradient so that each sub -pulse of the binomial 
pulse train is applied during portions of the oscillating 
gradient having the same polarity. The oscillating gradi- 
ent further includes at least one extra gradient switch 
added to the end of the oscillating gradient for providing 
a balanced gradient for inherently mot ion -compensated 
slice selection. 

The method can be applied for performing motion- 
compensated spectral spatial selective magnetic reso- 
nance angiography by using a binomial pulse train hav- 
ing sub-pulse envelopes which are optimized for provid- 
ing slice variable tip angle producing a spatially varying 
slice excitation profile for a predetermined range of 
blood flow velocities and slab widths and a predeter- 
mined magnetic field strength. The spatially varying 
slice excitation profile can be a non-linearly increasing 
slice excitation profile. The methods provides an im- 
proved dynamic range of the MR I images and improved 
small vessel visualization in angiography. 
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Description 

FIELD OF THE INVENTION 

The present, invention relates generally to the field s 
of medical imaging and more particularly relates to the 
field of magnetic resonance imaging (MRI). 

BACKGROUND OF THE INVENTION 

A unique property of magnetic resonance imaging 
(MRI) is the ability to selectively image different chemi- 
cal species by virtue of what is known as the chemical 
shift phenomenon. For example, in the human body the 
fat exhibits a resonance which is separated from the wa- 
ter resonance by approximately 3.5 ppm. At a field 
strength of 2 Tesla this equates to a frequency separa- 
tion of about 280 Hz in the NMR spectrum. Silicon ex- 
hibits a chemical shift of approximately 5 ppm. 

It is common that a "water only" image is required 
in which the fat or silicon resonance must be either sat- 
urated prior to each data acquisition or not excited at all. 

Suppression Techniques 

Techniques of fat suppression, generally referred to 
as FATS AT hereinafter, are well known in the prior art. 
Fat suppression is usually achieved by placing a narrow 
band spectral suppression pulse before the imaging se- 
quence. This pre-pulse is quickly followed by the imag- 
ing sequence so that the fat protons do not have time to 
relax back to their equilibrium magnetization which re- 
mains dispersed (saturated) and unable to contribute 
signal to the image. 

Reference is now made to Fig..1 which is a sche- 
matic graph Illustrating a prior art 3D gradient-echo mo- 
tion compensated sequence having a frequency selec- 
tive simple gaussian FATSAT spectral suppression 
pulse which is centered on the fat frequency. All the hor- 
izontal axes of Fig. 1 represent time. The double-head- 
ed arrow labeled Tp represents the presaturation se- 
quence and the double-headed arrow labeled TR rep- 
resents the imaging sequence. The graph labeled RF 
represents the RF imaging sequence which includes a 
gaussian spectral suppression pre-pulse 2 followed by 
an imaging RF pulse 3 having a flip angle a. 

The vertical axis (not shown) of the graph labeled 
RF represents the amplitude of the RF pulses. The 
graphs labeled Gs, Ge and Gv represent the slice se- 
lection gradient, the phase encoding gradient and the 
viewing gradient sequences, respectively, the vertical 
axes (not shown) of the graphs labeled Gs, Ge and Gv 
represent the gradient magnetic field amplitude. The 
slice selection gradient sequence includes a three- 
lobed motion compensated gradient 5, a phase encod- 
ing gradient 6 and a rewinder gradient 7. The gradient 
pulse referenced 4 is a spoiler pulse which is part of the 
presaturation sequence Tp. The phase encoding axis 



sequence includes a phase encoding gradient 8 and a 
rewinder gradient 9. The viewing gradient sequence (re- 
adout gradient) includes a readout gradient 10. 

The graph labeled S is the signal graph. The vertical 
axis (not shown) of the signal graph represents the RF 
echo signal amplitude. The signal graph includes the RF 
echo signal 11. 

Reference is now made to Fig. 2 which is a sche- 
matic graph representing a prior art NMR spectrum of 
fat and water protons on which a fat gaussian spectral 
suppression pulse is superimposed. The horizontal axis 
represents the chemical shift in parts per million (PPM) 
units. The curve labeled 14 represents the absorption 
spectrum of fat and water protons. The arrow labeled W 
indicates the peak absorption of the water protons at 0 
ppm and the arrow labeled F indicates the peak of ab- 
sorption of the fat protons which is shifted by 3.5 ppm 
relative to the peak absorption of the water protons. The 
excitation spectrum of a typical fat gaussian spectral 
suppression pulse 16 is superimposed on the absorp- 
tion spectrum curve 14. The gaussian suppression 
pulsel6 is centered at the peak F and will thus selectively 
excite the fat protons without substantial excitation of 
the water protons. Suppression pulses using the Sine 
function are also known in the art. 

Gaussian and Sine type suppression pulses are re- 
quired to be long in duration in order to achieve a suit- 
ably narrow spectral selection. At a field strength of 2 
Tesla a typical Sine suppression pulse may take up to 
26 ms. 

Reference is now made to Fig. 3 which is a sche- 
matic graph illustrating a prior art 3D gradient-echo im- 
aging sequence having a frequency selective binomial 
FATSAT spectral suppression presaturation pulse which 
is centered on the water frequency. 

The fat suppression sequence of Fig. 3 uses a 1 -3- 
3-2 binomial suppression pulse which is centered 
around the water frequency. 

The use of binomial pulse suppression techniques 
is disclosed in an article appearing in The Journal of 
Magnetic Resonance, entitled "Solvent Suppression in 
Fourier Transform Nuclear Magnetic Resonance" by P. 
J. Hore (Vol. 55, 1983, pp. 283 - 300) incorporated here- 
in by reference. 

It is noted that for a better understanding, like com- 
ponents are designated by like reference numerals 
throughout the various figures. 

The gradient-echo sequence of Fig. 3 is similar to 
the gradient -echo sequence of Fig. 1 , except that the 
RF imaging sequence which includes a gaussian spec- 
tral suppression pre-pulse 2 of Fig. 3 includes a 1 -3-3-1 
binomial suppression pulse 1 instead of the gaussian 
spectral suppression pre-pulse 2 of Fig. 1 . The 1 -3-3-2 
binomial suppression pulse 1 includes four sub-pulses 
1 A, 1 B, 1 C and 1 D which are separated from each other 
by a pulse separation interval t. Similar to the graphs of 
Fig. 1 , the horizontal axes of the graphs of Fig. 3 repre- 
sent time, the vertical axes (not shown) of the graphs 
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labeled, Gs, Ge, Gv represent gradient magnetic field 
amplitude, the vertical axis (not shown) of the graph la- 
beled RF represents the amplitude of the RF pulses and 
the vertical axis (not shown) of the graph labeled S rep- 
resents the RF echo signal amplitude. 

By choosing the appropriate pulse separation inter- 
val x (dependent upon field strength and chemical spe- 
cies) the binomial pulse 1 exhibits a null excitation at the 
water frequency which rises to a 90° excitation at the fat 
frequency. 

Reference is now made to Fig. 4 which is a sche- 
matic graph illustrating the theoretical excitation spec- 
trum of a prior art 1-3-3-1 binomial suppression pulse 
as a function of frequency offset from the transmitter fre- 
quency. The vertical axis of the graph represents the 
transverse magnetization Mxy wherein full scale corre- 
sponds to complete conversion of Z magnetization into 
XY magnetization. The horizontal axis represents the 
offset From the transmitter frequency. The transverse 
magnetization curve 40 has a flat excitation null 42 
around the transmitter frequency. . 

The binomial pulse sequence shown in Fig. 3 has 
a total duration of approximately 5.4 ms at 2 Tesla. This 
is somewhat shorter than the gaussian pulse 2 of Fig. 1 
but requires a high RF power because of the short "hard" 
pulses. 

Suppression techniques generally extend the min- 
imum TR that can be used and result in a reduction in 
the number of slices that can be imaged in a multi-slice 
sequence. They are also limited when short TR's are 
required since rapid, repeated, and incomplete, satura- 
tion of the fat frequency inevitably leads to a build up of 
coherent fat signal resulting in image artifacts. 

Spectral -spatial excitation methods 

Methods of Spectral-Spatial Excitation use a care- 
fully designed RF modulation in the presence of an os- 
cillating gradient. The result is a simultaneous selection 
along one spatial axis and the frequency spectrum. The 
use of Spectral-spatial excitation is disclosed in an arti- 
cle appearing in Magnetic Resonance in Medicine, en- 
titled "Simultaneous Spatial and Spectral selective Ex- 
citation" by Craig H. Meyer et al. (Vol. 15, 1990, pp. 287 
- 304), incorporated herein by reference. 

Fig. 5 is a schematic graph illustrating imaging se- 
quences designed for a prior art spectral -spatial excita- 
tion method. The horizontal axes of all the graphs of Fig. 
5 represent time. The graph labeled RF represents the 
RF "fat free" imaging sequence. The vertical axis (not 
shown) of the graph labeled RF represents the ampli- 
tude of the RF. The graphs labeled G z , G x and G y rep- 
resent the Gradient sequences wherein the horizontal 
axis represents time and the vertical axes(not shown) 
represent the magnetic field amplitude. The gradient se- 
quence labeled G z is a modulated slice selection gradi- 
ent. The gradient sequences labeled Gx and G y are spi- 
ral readout gradients. Each of the gradients G 2 , G x and 



G y is shown as having a rephasing pulse at the far end 
of the gradient pulse trains. 

The frequency of the modulated gradient is calcu- 
lated so that, when centered on the water resonance, 

s an excitation null occurs at the fat resonance. In this way 
only the water resonance is excited. This kind of pulse 
is usually incorporated directly into the imaging se- 
quence since it is designed to select the desired slice 
profile only at the water frequency. 

10 Reference is now made to Fig. 6 which is a sche- 
matic graph illustrating a prior art rapid gradient-echo 
pulse sequence using a spectral-spatial pulse with 
gaussian k-space varying along both k z and k^ as dis- 
closed by Meyer et al. The horizontal axes of all graphs 

is represent time. The graphs labeled RFI and RFQ rep- 
resent the real and the imaginary components of the RF 
"fat free" imaging sequence, respectively. In the partic- 
ular example of Fig. 6, there is no imaginary component 
as illustrated by the flat trace of the curve labeled RFQ. 

20 The vertical axis (not shown) of the graphs labeled RFI 
and RFQ represent RF signal amplitude. The graphs la- 
beled G z , G x and G y represent the Gradient sequences 
wherein the horizontal axis represents time and the ver- 
tical axes (not shown) represent the magnetic field am- 

25 plitude. The gradient sequence labeled G z is a modu- 
lated slice selection gradient. The gradient sequences 
labeled G x is a readout gradient. The gradient sequence 
G y is a phase encoding gradient. This pulse sequence 
results in compact spatial and spectral slice profiles (not 

so shown) which are gaussian in shape in the small-tip-an- 
gle regime. 

Spectral -spatial techniques have the advantage of 
exciting only the chemical species of interest. Because 
of this, no sacrifice is necessary on the repetition time 
35 (TR). 

However, spectral-spatial pulses are limited by gra- 
dient performance and are especially limited for low field 
applications where they are prohibitively long in dura- 
tion. Additionally, careful optimization is required to en- 
40 sure good spectral selection. 

SUMMARY OF THE INVENTION 

There is therefore provided, in accordance with an 
<5 embodiment of the present invention, a method for per- 
forming motion -compensated spectral-spatial selective 
magnetic resonance imaging. The method includes ap- 
plying a radio frequency (RF) spectral-spatial excitation 
sequence to a region of a body. The excitation sequence 
so includes a binomial slice-selective pulse train including 
at least two sub-pulses. The binomial pulse train also 
provides spectrally selective excitation of a selected res- 
onance frequency at a predetermined field strength. The 
method further includes applying to the same region of 
55 the body an oscillating slice selection gradient so that 
each sub-pulse of the binomial pulse train is applied dur- 
ing portions of the oscillating gradient having the same 
polarity. The oscillating gradient further includes at least 
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one extra gradient switch added to the end of the oscil- 
lating gradient for providing inherently motion compen- 
sated slice selection. The method thereby provides an 
improved dynamic range of MRI imaging. 

Further, in accordance with another embodiment of 
the present invention, the selected resonance frequen- 
cy is the water resonance frequency for the predeter- 
mined field strength thereby reducing artifacts due to ex- 
citation of the fat resonance frequency and providing 
substantially fat free images. 

Further still, in accordance with yet another embod- 
iment of the present invention, the binomial pulse train 
of the RF excitation sequence is phase modulated in or- 
der to shorten the duration of the slice selection process 
and minimize the TE duration. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the binomial pulse train 
is a 1-2-1 binomial pulse train. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the binomial pulse train 
is a 1-3-3-1 binomial pulse train. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the envelope of each sub- 
pulse of the RF excitation sequence is modified to cor- 
rect for the effect of the time varying of the waveform of 
the oscillating slice selection gradient, by using non-lin- 
ear sampling methods or other suitable methods of 
pulse envelope correcting, thereby maintaining the fidel- 
ity of the slice excitation profile. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the spectrally selective 
excitation at a selected resonance frequency is pro- 
duced by selecting the spacing interval between the 
sub-pulses of the binomial pulse train in accordance 
with the predetermined field strength. 

Furthermore, in accordance with another embodi- 
ment of the present invention, each sub-pulse of the bi- 
nomial pulse train is applied during portions of the os- 
cillating gradient having identical gradients. 

Furthermore, in accordance with another embodi- 
ment of the present invention, each sub-pulse of the bi- 
nomial pulse train is applied during positive gradient por- 
tions of the oscillating gradient. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the selected resonance 
frequency is the water resonance frequency at the pre- 
determined field strength. 

There is also provided, in accordance with another 
embodiment of the present invention, a method for per- 
forming motion -compensated spectral-spatial selective 
magnetic resonance angiography. The method includes 
applying a radio frequency (RF) spectral-spatial excita- 
tion sequence to a region of a body. The excitation se- 
quence includes a binomial pulse train having at least 
two sub-pulses. Each of the sub-pulses of the binomial 
pulse train has an envelope optimized for providing 
slice-variable tip angle producing a spatially varying 
slice-excitation-profile for a predetermined range of 



blood flow velocities and slab widths and a predeter- 
mined magnetic field strength. The binomial pulse train 
also provides spectrally selective excitation of a select- 
ed resonance frequency at the predetermined field 

5 strength. The method also includes applying to the 
same region of the body an oscillating slice selection 
gradient so that each sub-pulse of the binomial pulse 
train is applied during portions of the oscillating gradient 
having the same polarity. The oscillating gradient further 

10 includes at least one extra gradient switch added to the 
end of the oscillating gradient for providing inherently 
flow compensated slice selection. The method provides 
an improved dynamic range of the MR angiography and 
an improved small vessel visualization. 

is Furthermore, in accordance with another embodi- 
ment of the present invention, the spatially varying slice 
excitation profile is a non-linearly increasing slice exci- 
tation profile. 

Further yet, in accordance with another embodi- 

20 ment of the present invention, the selected resonance 
frequency is the water resonance frequency for the pre- 
determined field strength thereby reducing artifacts due 
to excitation of the fat resonance frequency and provid- 
ing substantially fat free images. 

25 Furthermore, in accordance with another embodi- 
ment of the present invention, the binomial pulse train 
of the RF excitation sequence is phase modulated in or- 
der to shorten the duration of the slice selection process 
and minimize the TE duration. 

30 Furthermore, in accordance with another embodi- 
ment of the present invention, the binomial pulse train 
is a 1-2-1 binomial pulse train. 

Furthermore, in accordance with yet another em- 
bodiment of the present invention, the binomial pulse 

35 train is a 1 -3-3-1 binomial pulse train. 

Furthermore, in accordance with another embodi- 
ment of the present invention, the envelope of each sub- 
pulse of the RF excitation sequence is modified to cor- 
rect for the effect of the time varying of the waveform of 

40 the oscillating slice selection gradient, thereby maintain- 
ing the fidelity of the slice excitation profile. 

Furthermore, in accordance with another embodi- 
. ment of the present invention, the spectrally selective 
excitation at a selected resonance frequency is pro- 

45 duced by selecting the spacing interval between the 
sub-pulses of the binomial pulse train in accordance 
with the predetermined field strength. 

Furthermore, in accordance with another embodi- 
ment of the present invention, each sub-pulse of the bi- 

50 nomial pulse train is applied during portions of the os- 
cillating gradient having identical gradients. 

Furthermore, in accordance with another embodi- 
ment of the present invention, each individual pulse of 
the binomial pulse train is applied during positive gradi- 

55 ent portions of the oscillating gradient. 

Finally, in accordance with another embodiment of 
the present invention, the selected resonance frequen- 
cy is the water resonance frequency at the predeter- 
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mined field strength. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of exam- s 
pie only, with reference to the accompanying drawings, 
wherein: 



1 -2-1 BINSS pulse sequence designed for fat-free 
angiography using motion compensated binomial- 
spectral-spatial-slice-variable-tip-angle RF slice 
selection in accordance with another embodiment 
of the present invention; 

Fig. 16 is a schematic three dimensional graph il- 
lustrating the spectral and spatial selectivity of the 
1-2-1 BINSS pulse of Fig. 15; 
Fig. 17 is a schematic 2-dimensiona! contour graph 
illustrating the spectral and spatial selectivity of the 
1-2-1 BINSS pulse of Fig. 15; 
Fig. 18 is a schematic graph illustrating the simulat- 
ed slice profile curve generated by the 1 -2-1 BINSS 
pulse of Fig. 15; 

Fig. 19 is a schematic graph illustrating the simulat- 
ed spectral response curve of the 1-2-1 BINSS 
pulse of Fig. 15; and 

Fig. 20 is a schematic graph illustrating a full imag- 
ing sequence for fat-free angiography using a mo- 
tion compensated 1-2-1 binomial-spectral-spatial 
-slice-variable-tip-angle RF slice selection in ac- 
cordance with another embodiment of the present 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 

While the invention has been described with re- 
spect to a limited number of embodiments, it will be ap- 
preciated that many variations, modifications and other 
applications of the invention may be made. 

Binomial Spectral Spatial Excitation 

The binomial spectral-spatial excitation method, re- 
ferred to as the BINSS method hereinafter, combines 
the simplicity and compact nature of binomial pulses 
with the principle of spectral -spatial excitation to pro- 
duce a slice selection pulse which has excellent selec- 
tivity in both the spectral and spatial domains. 

Reference is now made to Fig. 7 which is a sche- 
matic graph illustrating a 1 -2-1 BINSS RF pulse together 
with the slice selection gradient selected for that pulse, 
in accordance with an embodiment of the present inven- 
tion. In the graph labeled RF of Fig. 7 the curve labeled 
50 represents the normalized magnitude of the binomial 
RF fat free imaging pulse. The horizontal axis repre- 
sents time and the vertical axis represents the normal- 
ized amplitude of the RF sequence. In the graph of Fig. 
7 labeled Gs, the curve labeled 54 represents the mod- 
ulated slice selection gradient sequence. The horizontal 
axis represents time and the vertical axis represents the 
normalized amplitude of the slice selection gradient. 
The phase encoding gradient and the viewing gradient 
are not shown. 

The correct sub-pulse spacing x for "water only" im- 
aging at a selected field strength is given by: 



Fig. 1 is a schematic graph illustrating a prior art 3D 
gradient-echo motion compensated sequence hav- to 
ing a frequency selective simple gaussian FATSAT 
spectral suppression pulse; 
Fig. 2 is a schematic graph representing a prior art 
NMR spectrum of fat and water protons on which a 
fat gaussian spectral suppression pulse is 15 
superimposed ; 

Fig. 3 is a schematic graph illustrating a prior art 3D 
gradient-echo imaging sequence having a frequen- 
cy selective binomial FATSAT spectral suppression . 
presaturation pulse centered on the water frequen- 20 
cy; 

Fig. 4 is a schematic graph illustrating the theoreti- 
cal excitation spectrum of a prior art 1 -3-3-2 Dino ~ 
mial suppression pulse as a function of frequency 
offset from the transmitter frequency; 2s 
Fig. 5 is a schematic graph illustrating imaging se- 
quences designed for a prior art spectral-spatial ex- 
citation method; 

Fig. 6 is a schematic graph illustrating a prior art 
rapid gradient-echo pulse sequence using a spec- 30 
tral-spatial pulse with gaussian k-space varying 
along both k z and k^ ; 

Fig. 7 is a schematic graph illustrating a 1 -2-1 bino- 
mial spectral-spatial (BINSS) RF pulse together 
with the slice selection gradient selected for that 35 
pulse; 

Fig. 8 is a schematic graph illustrating the slice pro- 
file of the BINSS pulse of Fig. 7; 
Fig. 9 is a schematic graph illustrating the spectral 
response of the BINSS pulse of Fig. 7; *o 
Fig. 10 is a schematic graph illustrating a motion 
compensated 1-3-3-1 BINSS pulse and a slice se- 
lection gradient designed for fat free imaging in ac- 
cordance with an embodiment of the present 
invention ; 45 
Fig. 11 is a schematic three dimensional graph il- 
lustrating the spectral and spatial selectivity of the 
RF pulse sequence of Fig. 10 ; 
Fig. 12 is a schematic 2-dimensional contour graph 
illustrating the spectral and spatial selectivity of the so 
RF pulse sequence of Fig. 10; 
Fig. 1 3 is a schematic graph illustrating the simulat- 
ed slice profile generated by the 1-3-3-1 BINSS 
pulse of Fig. 10; 

Fig. 14 is a schematic graph illustrating the spectral ss 
response curve of the 1-3-3-1 BINSS pulse of Fig. 
10; 

Fig. 15 is a schematic graph illustrating a complex 
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T 2. AW 

wherein Af FAT is the frequency separation be- 
tween the fat and the water resonance for the selected 
field strength. 

The calculated value of x, at a field strength of 2 
Tesla, is 1784 micro-seconds. For the binomial 1-2-1 RF 
sequence of Fig. 7, three sub-pulses are required lead- 
ing to a total pulse duration of 5.352 ms as opposed to 
the 1 0.7 ms required for prior art spectral-spatial pulses. 

At a field strength of 0.5 Tesla (not shown) the pulse 
spacing is 6756 micro-seconds. For a 1-2-1 binomial 
pulse at a field strength of 0.5 Tesla the total pulse du- 
ration is 20 ms, in contrast to the 40.5 ms required for 
prior art spectral-spatial pulses. 

Non-Linear Correction of Sub-Pulses 

Because the RF pulses are applied under a time 
varying gradient waveform, it is necessary to correct the 
sub-pulse shape to take the gradient waveform into ac- 
count so that fidelity in the slice profile is maintained. 

Binomial spectral-spatial pulse sequences are 
shorter in duration thus enabling shorter TE's in imag- 
ing. Additionally, a saving in gradient performance by a 
factor of 2, for the same slice thickness, is obtained 
when using the BINSS pulse over prior art. Furthermore, 
the Bl NSS pulse can be made arbitrarily longer or short- 
er by using other binomial series. For example, a 1 -3-3-1 
binomial pulse train can be used at high magnetic field 
strengths such as 2 tesla, while a 1 : 2-1 binomial pulse 
train can be used for spectral-spatial excitation at lower, 
field strengths, such as 0.5 Tesla. 

Spatial and Spectral Response of BINSS pulse 

Reference is now made to Figs. 8 and 9. Which are 
schematic graphs illustrating the slice profile and the 
spectral response of the BINNS pulse of Fig. 7. 

Fig. 8 illustrates the simulated slice profile curve 56 
generated by the 5.532 ms 1 -2-1 BINSS pulse of Fig. 7. 
The horizontal axis represents the distance along the 
spatial selection coordinate and the vertical axis repre- 
sents the normalized excitation amplitude. It is noted 
that the slice profile is symmetrical. 

Fig. 9 illustrates the spectral response curve 58 of 
the BINSS pulse of Fig. 7. The horizontal axis repre- 
sents the frequency in Hz relative to the water reso- 
nance frequency (at a field strength of 2 Tesla) and the 
vertical axis represents the normalized excitation ampli- 
tude. The peak excitation is at the Water resonance fre- 
quency. 

Motion Compensated spectral -spatial BINSS pulse 

Reference is now made to Fig. 1 0 which is a sche- 



matic graph illustrating a motion compensated spectral- 
spatial 1-3-3-1 BINSS pulse and a slice selection gradi- 
ent designed for fat free imaging in accordance with an 
embodiment of the present invention. The graph labeled 

5. RF illustrates the RF BINSS pulse and the graph labeled 
Gs illustrates the slice selection gradient. The curve la- 
beled 55 represents the real component of the binomial 
RF fat free imaging sequence. It is noted that, the pulse 
sequence has a null imaginary component (not shown). 

10 The horizontal axis represents time and the vertical axis 
represents the normalized amplitude of the RF se- 
quence for the graph labeled RF. The curve labeled 57 
represents the modulated slice selection gradient se- 
quence. In the graph labeled Gs, the horizontal axis rep- 

15 resents time and the vertical axis represents the normal- 
ized gradient amplitude along the Z axis. The gradient 
is switched at double the rate of the corresponding gra- 
dient of Fig. 7 and the sub-pulses of the binomial RF 
pulse sequence are applied during the positive gradient 

20 switches only. 

The phase encoding and the viewing gradients for 
the RF pulse of Fig. 10 are not shown. 

The end of the gradient waveform 57 is balanced 
using the prior art method of gradient moment nulling 

25 (GMN), so that motion compensation is achieved. Since 
the RF sub-pulses are applied during the same gradient 
polarity each time there is no need to modulate the 
phase of the pulses in order to achieve shifted slices. 
Thus, in contrast to the 1 -2-1 BINSS pulse of Fig. 7 hav- 

30 ing sub-pulses applied during opposite gradient polari- 
ties, the sub-pulses of the 1-3-3-1 BINSS pulse of Fig. 
1 0 are applied during the same gradient polarities, thus, 
taking advantage of the inherent motion -compensation 
of the switched gradient . 

35 it is noted that the sub-pulses of the binomial RF 
pulse sequence can also be applied during the negative 
gradient switches as long as they are applied during the 
same gradient polarity each time. 

Reference is now made to Figs 11 and 12. Fig. 11 

40 is a schematic three dimensional graph illustrating the 
spectral and spatial selectivity of the RF pulse sequence 
of Fig. 10. The horizontal axes labeled X, and Y repre- 
sent the relative frequency (spectral axis) and the slice 
excitation profile (spatial axis), respectively. The vertical 

45 axis represents the final magnetization Mxy of the 7. 1 36 
ms 1-3-3-1 BINSS pulse for water-only excitation at 2 
Tesla. 

Fig. 12 is a schematic 2-dimensional contour graph 
illustrating the spectral and spatial selectivity of the RF 
50 pulse sequence of Fig. 1 0. The horizontal axes labeled 
X, and Y represent the relative frequency (spectral axis) 
and the slice excitation profile (spatial axis), while the 
contour lines represent lines of equal final magnetiza- 
tion Mxy. 

55 Reference is now made to Figs. 13 and 14. Fig. 13 
is a schematic graph illustrating the simulated slice pro- 
file curve referenced 60 generated by the 7.136 ms 
1 -3-3-1 BINSS pulse of Fig. 10. The horizontal axis rep- 



6 



11 



EP 0 871 039 A2 



12 



resents the distance along the spatial selection axis and 
the vertical axis represents the normalized excitation 
amplitude. It is noted that the slice profile is symmetrical. 

Fig. 14 is a schematic graph illustrating the spectral 
response curve referenced 62 of the BINSS pulse of Fig. 
10. The horizontal axis represents the frequency in Hz 
relative to the water resonance frequency (at a field 
strength of 2 Tesla) and the vertical axis represents the 
normalized excitation amplitude. The peak excitation is 
at the Water resonance frequency. 

Motion Compensated slice variable tip angle BINSS 
pulse angiography 

Reference in now made to Fig. 15 which is a sche- 
matic graph illustrating a complex 1-2-1 BINSS pulse 
sequence designed for fat-free angiography using mo- 
tion compensated binomial-spectral : spatial-slice-varia- 
ble-tip-angle RF slice selection in accordance with an- 
other embodiment of the present invention. The curve 
64 of the graph labeled RFR represents the real com- 
ponent of the binomial RF fat free imaging sequence. 
The horizontal axis represents time and the vertical axis 
represents the normalized amplitude of the real compo- 
nent of the complex RF pulse. 

The curve 66 of the graph labeled RFl represents 
the imaginary component of the binomial RF fat free im- 
aging sequence. The horizontal axis represents time 
and the vertical axis represents the amplitude of the im- 
aginary component of the complex RF pulse. It is noted 
that, the curves 64 and 66 of the graphs labeled RFR 
and RFl, respectively, together represent the normal- 
ized complex RF pulse. 

The curve 68 of the graph labeled GS represents a 
balanced modulated slice selection gradient sequence 
wherein the horizontal axis represents time and the ver- 
tical axis represents the normalized gradient amplitude 
along the Z axis. It is noted that, while the slice selection 
gradient of the exemplary sequence of Fig. 1 5 is applied 
to the Z axis, the slice selection gradient can be applied 
to another spatial axis (not shown) such as the X or Y 
axis depending on the application. Additionally, more 
than one slice selection gradient can be applied to se- 
lected axis combinations (not shown) depending on the 
application. The phase encoding and the viewing gradi- 
ents for the RF pulse of Fig. 1 5 are not shown. 

Each of the three RF sub-pulses are complex puls- 
es designed to generated a slice profile in which the tip 
angle varies across the slice. The use of RF pulses for 
generating variable tip-angle slice profiles (referred to 
as SVTA pulse, hereinafter) is already widely used in 
inflow angiography. For example, the use of an SVTA 
pulse is disclosed in an article appearing in SMRM ab- 
stracts, entitled "Optimized Slab Profile for 3D-TOF An- 
giography" by Yuval Zur. (Abstract 960, 1 994), incorpo- 
rated herein by reference. 

Incorporating an SVTA pulse envelope in a 1-2-1 
BINSS pulse combined with a motion compensated 



slice selection gradient waveform provides an improved 
method of slice selection for performing fat free angiog- 
raphy. The optimized sub-pulses of the 1-2-1 RF pulse 
sequence of Fig. 15 were calculated using the method 

s disclosed by Yuval Zur. The advantages of using water 
only excitation in angiography are the greater dynamic 
range and improved small vessel visualization. 

Reference is now made to Figs 16 and 17. Fig. 16 
is a schematic three dimensional graph illustrating the 

io spectral and spatial selectivity of the RF pulse sequence 
of Fig. 15. The horizontal axes labeled X, and Y repre- 
sent the relative frequency (spectral axis) and the slice 
excitation profile (spatial axis), respectively. The vertical 
axis Z represents the final magnetization Mxy of the 

J5 5.352 ms 1 -2-1 BINSS pulse for water only excitation at 
a magnetic field strength of 2 Tesla. 

Fig. 17 is a schematic 2-dimensional contour graph 
illustrating the spectral and spatial selectivity of the RF 
pulse sequence of Fig. 15. The horizontal axes labeled 

20 X, and Y represent the relative frequency (spectral axis) 
and the slice excitation profile (spatial axis), while the 
contour lines represent lines of equal final magnetiza- 
tion Mxy. 

Reference is now made to Figs. 18 and 19. Fig. 18 

25 is a schematic graph illustrating the simulated slice pro- 
file curve referenced 70 generated by the 1-2-1 BINSS 
pulse of Fig. 15. The horizontal axis represents the dis- 
tance along the spatial selection axis and the vertical 
axis represents the normalized excitation amplitude. It 

30 is noted that the slice profile is non-symmetrical and in- 
creases nbn-linearly in the blood flow direction. The ar- 
row labeled FD represents the blood flow direction. 

Fig. 19 is a schematic graph illustrating the simulat- 
ed spectral response curve referenced 72 of the BINSS 

35 pulse of Fig. 15. The horizontal axis represents the fre- 
quency in Hz relative to the water resonance frequency 
(at a field strength of 2 Tesla) and the vertical axis rep- 
resents the normalized excitation amplitude. The peak 
of the spectral response curve 72 is at the water reso- 

40 nance frequency. It is noted that the excitation amplitude 
at the fat frequency which is offset from the water reso- 
nance by -280 Hz (at a field strength of 2 Tesla) is almost 
null indicating the selective excitation of water protons. 
Reference is now made to Fig. 20 which is a sche- 

45 matic graph illustrating a full imaging sequence for fat- 
free angiography using motion compensated 1 -2-1 bi- 
nomial-spectral-spatial-slice-variable-tip-angle RF slice 
selection in accordance with yet another embodiment of 
the present invention. All the horizontal axes of Fig. 20 

50 represent time. The double-headed arrow labeled TR 
represents the imaging sequence. The graph labeled 
RF represents the RF imaging sequence which includes 
a 1-2-1 binomial pulse train 80. The pulse train 80 in- 
cludes three sub-pulses 80A, 80B and 80C, having flip 

55 angles a, 2a and a, respectively. The sub-pulse sepa- 
ration interval is x. The vertical axis (not shown) of the 
graph labeled RF represents the amplitude of the RF 
pulses. 
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It is noted that, the RF pulse 80 can be a complex 
or a non-complex RF pulse. 

The graphs labeled Gs, Ge and Gv represent the . 
slice selection gradient, the phase encoding gradient 
and the viewing gradient sequences, respectively, the 
vertical axes (not shown) of the graphs labeled Gs, Ge 
and Gv represent the gradient magnetic field amplitude. 
The slice selection gradient sequence includes a bal- 
anced seven-lobed motion compensated gradient 82, a 
phase encoding gradient 86 and a rewinder gradient 88. 
The phase encoding axis sequence includes a phase 
selection gradient 90 and a rewinder gradient 92. The 
viewing gradient sequence (readout gradient) includes 
a readout gradient 94. 

The graph labeled S is the signal graph. The vertical 
axis (not shown) of the signal graph represents the RF 
echo signal amplitude. The signal graph includes the RF 
echo signal 96. 

It is noted that, while the exemplary embodiment of 
Fig. 20 uses a 1-2-1 binomial RF pulse, other suitable 
types of binomial RF pulses such as 1-3-3-1 or 1-3-3-1. 
can be used by suitably optimizing the gradient se- 
quences to conserve their motion compensating design. 

It will be appreciated by those skilled in the art that, 
while the embodiments disclosed hereinabove are ap- 
plied for the non-limiting example of spectrally selective 
water protons excitation ( with fat non -excitation), the 
method of the present invention can be applied to the 
selective excitation of other species such as fat protons 
(with water non -excitation) and is generally applicable 
for selective excitation of a desired species in the pres- 
ence of other interfering species having a resonance fre- 
quency which is close but not identical to the resonance 
frequency of the desired species. 

It will be appreciated by those skilled in the art, that 
the present invention is not limited to what have been 
particularly shown and described hereinabove. Rather, 
the scope of the present invention is defined by the 
claims which follow. 



Claims 

1. A method for performing motion -compensated 
spectral-spatial selective magnetic resonance im- 
aging (MRl), the method comprising: 

applying a radio frequency (RF) spectral-spa- 
tial excitation sequence to a region of a body, 
said excitation sequence comprising a binomial 
slice-selective pulse train, said binomial pulse 
train being comprised of at least two sub-puls- 
es, said binomial pulse train also providing 
spectrally selective excitation of a selected res- 
onance frequency at a predetermined field 
strength; 

applying to said region an oscillating slice se- 
lection gradient so that each sub-pulse of said 



binomial pulse train is applied during portions 
• of said oscillating gradient having the same po- 
larity, said oscillating gradient further compris- 
ing at least one extra gradient switch added to 
5 the end of said oscillating gradient for providing 

inherently motion compensated slice selection, 
for providing an improved dynamic range of 
said MRl imaging. 

10 2. The method of claim 1 wherein said selected reso- 
nance frequency is the water resonance frequency 
for said predetermined field strength thereby reduc- 
ing artifacts due to excitation of the fat resonance 
frequency and providing substantially fat free imag- 
es es. 

3. The method of claim 1 or claim 2 wherein said bi- 
nomial pulse train of said RF excitation sequence 
is phase modulated in order to shorten the duration 

20 of the slice selection process and minimize the TE 
duration. 

4. The method of any of claims 1 -3, wherein said bi- 
nomial pulse train is a 1-2-1 binomial pulse train. 

25 

5. The method of any of claims 1 -3, wherein said bi- 
nomial pulse train is a 1-3-3-1 binomial pulse train. 

6. The method of any of claims 1 -5, wherein the enve- 
30 lope of each sub-pulse of said RF excitation se- 
quence is modified to correct for the effect of the 
time varying of the waveform of said oscillating slice 
selection gradient, for maintaining the fidelity of said 
slice excitation profile. 

35 

7. The method of any of claims 1 -6 wherein said spec- 
trally selective excitation at a selected resonance 
frequency is produced by selecting the spacing in- 
terval between said sub-pulses of said binomial 

40 pulse train in accordance with said predetermined 
field strength. 

8. The method of any claims 1-7 wherein each sub- 
pulse of said binomial pulse train is applied during 

45 portions of said oscillating gradient having identical 
gradients. 

9. The method of claim 8 wherein each sub-pulse of 
said binomial pulse train is applied during positive 

so gradient portions of said oscillating gradient. 

10. The method of any of claims 1-9 wherein said se- 
lected resonance frequency is the water resonance 
frequency at said predetermined field strength. 

55 

11. A method for performing mot ion -compensated 
spectral-spatial selective magnetic resonance ang- 
iography, the method comprising: 
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applying a radio frequency (RF) spectral-spa- 
tial excitation sequence to a region of a body, 
said excitation sequence comprising a binomial 
pulse train, said binomial pulse train being com- 
prised of at least two sub-pulses, each of said s 
sub-pulses of said binomial pulse train having 
an envelope optimized for providing slice-vari- 19. 
able tip angle producing a spatially varying 
slice-excitation-profile for a predetermined 
range of blood flow velocities and slab widths 10 
and a predetermined magnetic field strength, 
said binomial pulse train also providing spec- 20. 
trally selective excitation of a selected reso- 
nance frequency at said predetermined field 
strength; 1S 
applying to said region an oscillating slice se- 
lection gradient so that each sub-pulse of said 21. 
binomial pulse train is applied during portions 
of said oscillating gradient having the same po- 
larity, said oscillating gradient further compris- 20 
ing at least one extra gradient switch added to 
the end of said oscillating gradient for providing 
inherently flow compensated slice selection, 
for providing an improved dynamic range of 
said MR angiography and improved small ves- 25 
sel visualization. 

12. The method of claim 11 wherein said spatially var- 
ying slice excitation profile comprises a non-linearly 
increasing slice excitation profile. 30 

13. The method of claim 11 or claim 12 wherein said 
selected resonance frequency is the water reso- 
nance frequency for said predetermined field 
strength thereby reducing artifacts due to excitation 35 
of the fat resonance frequency and providing sub- 
stantially fat free images. 

14. The method of any of claims 11-13 wherein said bi- 
nomial pulse train of said RF excitation sequence 40 
is phase modulated in order to shorten the duration 

of the slice selection process and minimize the TE 
duration. 

15. The method of any of claims 1 1 -1 4 wherein said bi- *s 
nomial pulse train is a 1-2-1 binomial pulse train. 

16. The method of any of claims 11-14 wherein said bi- 
nomial pulse train is a 1-3-3-1 binomial pulse train. 

so 

17. The method of any of claims 11-16 wherein the en- 
velope of each sub-pulse of said RF excitation se- 
quence is modified to correct for the effect of the 
time varying of the waveform of said oscillating slice 
selection gradient, for maintaining the fidelity of said 55 
slice excitation profile. 

18. The method of any of claims 11-17 wherein said 



spectrally, selective excitation at a selected reso- 
nance frequency is produced by selecting the spac- 
ing interval between said sub-pulses of said bino- 
mial pulse train in accordance with said predeter- 
mined field strength. 

The method of any of claims 11-18 wherein each 
sub-pulse of said binomial pulse train is applied dur- 
ing portions of said oscillating gradient having iden- 
tical gradients. 

The method of claims 19 wherein each individual 
pulse of said binomial pulse train is applied during 
positive gradient portions of said oscillating gradi- 
ent. 

The method of any of claims 11-20 wherein said se- 
lected resonance frequency is the water resonance 
frequency at said predetermined field strength. 
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(57) A method for performing motion-compensated 
spectral-spatial selective magnetic resonance imaging. 
The method includes applying a radio frequency (RF) 
spectral-spatial excitation pulse sequence to a region of 
a body. The excitation sequence includes a binomial 
pulse train including at least two sub-pulses. The bino- 
mial pulse train also provides spectrally selective exci- 
tation of a selected resonance frequency at the prede- 
termined field strength. The method also includes ap- 
plying to the same body region an oscillating slice se- 
lection gradient so that each sub-pulse of the binomial 
pulse train is applied during portions of the oscillating 
gradient having the same polarity. The oscillating gradi- 
ent further includes at least one extra gradient switch 
added to the end of the oscillating gradient for providing 
a balanced gradient for inherently motion -compensated 
slice selection. 

The method can be applied for performing motion- 
compensated spectral spatial selective magnetic reso- 
nance angiography by using a binomial pulse train hav- 
ing sub-pulse envelopes which are optimized for provid- 
ing slice variable tip angle producing a spatially varying 
slice excitation profile for a predetermined range of 
blood flow velocities and slab widths and a predeter- 
mined magnetic field strength. The spatially varying 
slice excitation profile can be a non-linearly increasing 
slice excitation profile. The methods provides an im- 
proved dynamic range of the MR I images and improved 
small vessel visualization in angiography. 
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